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Analytical solutions for geodesics in black hole 
spacetimes 


Eva Hackmann and Claus Lammerzahl 


Abstract We review the analytical solution methods for the geodesic equations in 
Kerr-Newman-Taub-NUT-de Sitter spacetimes and its subclasses in terms of elliptic 
and hyperelliptic functions. A short guide to corresponding literature for general 
timelike and lightlike motion is also presented. 


1 Introduction 


Black holes belong to the most fascinating objects in astrophysics and are well suited 
to explore the regime of strong gravity. We consider here black holes which are 
described by the six parameter family of Kerr-Newman-Taub-NUT-de Sitter space- 
times including mass, rotation, electric and magnetic charge, gravitomagnetic mass 
(or NUT charge), and the cosmological constant. Maybe the best way to explore the 
gravitational field of such objects is through the observation of the motion of small 
massive particles and light, which is described by the geodesic equation. The com¬ 
plete set of orbits can best be analyzed using analytical methods. Already in 1931, 
Hagihara ESI used Weierstrass elliptic functions to analytically solve the geodesic 
equation in Schwarzschild spacetime. Later, Darwin lHISl solved the geodesic equa¬ 
tions in Schwarzschild spacetime in terms of Jacobian elliptic functions. These 
methods and their generalization to hyperelliptic functions can be used to solve the 
geodesic equation in the six parameter spacetime under consideration. Although this 
requires only mathematics from the 19th century, surprisingly the geodesic equa¬ 
tions in Schwarzschild-de Sitter spacetime were analytically solved only in 2008 
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muni. Here we will review these solution methods and provide a short literature 
guide. 


2 The Kerr-Newman-Taub-NUT-de Sitter space-time 


We consider here spherically or axially symmetric spacetimes with up to six param¬ 
eters, which are part of the family of type D Plebariski-Demiahski spacetimes. The 
metric is given by 0|23l 


n2 

P J..2 


ds^ =—^ (dt — {a sin^ 0+2n COS 9)d(p) - dr' 


^0 


P jq2 


- j sin"^ 9{adt — ir^ + n‘^)d(pY - d9'^ (1) 


with 

= r^+ (« —acos0)^ , 

Ag = I + ^a^Acos^ 9 — jAancos9 , (2) 

Ar = r^ — 2Mr + + qI + + (6«^ + + 3(fl^ — , 

where M is the mass, a=J/M the specific angular momentum, A the cosmological 
constant, n is the gravitomagnetic mass, Qe is the electric, and the magnetic 
charge of a gravitating source. 


3 Analytical solution methods 

The motion of test particles in the spacetime metric O is given by the geodesic 
equation 


dh^' 

ds^ 


+r^ 


dx'’ dxP 
ds ds 


( 3 ) 


where (F^vp) are the Christoffel symbols and s is an affine parameter. As the metric 
([1]) is axially symmetric there exists the two constants of motion 


E = = (4) 

connected to the Killing vectors and which can be interpreted as the en¬ 
ergy and the specific angular momentum in direction of the symmetry axes. Here u 
denotes the four-velocity. If we have even spherical symmetry, i.e. for a = n = 0, 
these two constants of motion together with the restriction to the equatorial plane. 
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which is without loss of generality, and the normalization g^ytix^dx'' = e, where 
e = 0 for light and £ = 1 for massive particles, is sufficient to separate the geodesic 
equation. This yields then the differential equation 



2 




( 5 ) 


and the energy and angular momentum take the form 


E =gn 


dt 

ds 


L = L- = k 


ds 


(6) 


If the spacetime is axially symmetric, we cannot in general restrict the motion of 
test particles to the equatorial plane. Therefore, we need an additional constant of 
motion. In 1968 Carter a surprisingly found such a constant of motion, which can 
be derived as a separation constant and is not connected to an obvious symmetry of 
the spacetime. With this constant the geodesic equation can again be separated and 
we get equations of motions in the form 


P 


4 



P 


4 



((r^ + + n^)E — aL,)^ 

Aq{K — e{n — a cos 6)^) 


— A,-{Er^ +-^)! 

(£’(asin^0 + 2ncos0) 
sin^ 0 




( 7 ) 

( 8 ) 


where K is the Carter constant. If n = 0 then K = {aE — corresponds to mo¬ 
tion restricted to the equatorial plane. From the expression of energy and angular 
momentum (IHi we get the additional equations 


2 d(p a ,, 2 2 t \ £'(2ucos0+asin^0)-L- 

p^^ = —{{r^ + a^ + n^)E-aL,)-^ - . . (9) 

ds 


Ag sin 0 


2 “' + + 2 , 2 


,dt 


^ ds 


Ar 


+ n^)E — aL^) 


asirr 9 + 2ncos6 
Ag sin^ 0 


(£’(2ncos 0 +asin^ 0) —L^). 


( 10 ) 


Note that the equations (|7]i and ® are still coupled by the factor p^. This issue was 
solved by Mino in 2003 ll25l by introducing a new affine parameter A defined by 

ds _ ,^2 

dX ~ P ■ 

If we now consider in the spherically symmetric case the differential equations 
for r{(p). 



2 






(11) 
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we see that on the right hand side we have a polynomial of degree three or four, if 
the cosmological constant vanishes, and of degree five or six in general. The same 
holds for the differential equations 0 and (| 8 ]l. This kind of differential equations 
can be solved in terms of elliptic functions if the polynomial is of degree three or 
four, and in terms of hyperelliptic functions if it is of degree five or six. 


3.1 Solutions in terms of elliptic functions 


For differential equations of the general type 

2 


(^) =-^ 3 , 4 ( x ), xiyo) 


= xo 


( 12 ) 


where F 3 4 is a polynomial of degree three or four, there are basically two (equiva¬ 
lent) solution methods based on the Jacobian elliptic function sn and on the Weier- 
strass elliptic function p. The first can be defined as the inverse of an elliptic inte¬ 
gral. 


z = 


rw 

Jo 


dt 


x/(l -f2)(l -k2f2) 


where 0 < k < 1 is the modulus, w G [ 0 , 1 ], 
function is given as a series 


P(z;2®i,2®2) = 


1 




^ sn(z;k) = w (13) 

and z £ R. The Weierstrass elliptic 


_LV 

(z (OnmY oflm J 


(14) 


where 2 tUi, 2 a >2 are the periods of p (^ ^ R) and Q = {(0„m G C\(Onm = 2n(Oi -f 
2ma)2}. It solves the initial value problem (see e.g. ll24l[T9l ) 

/ 7 \ 2 

=4x^-.g2X-g3, x(0)=°o (15) 

where g 2 = > ^3 = ^^0'Lco„„,en Note that both sn and pcan be 

written in terms of the Riemann theta function 


0 [tv -I-w](z; t) = ^ exp(;r/(m4-v)'(T(m-|-v)-|-2z + 2w)), (16) 

mezs 


where z G C^, T is a g x g symmetric matrix with positive definite imaginary part, 
Tv + w G C® is the characteristic, and g is the genus, here g = 1. 

The general differential equation (fT^ can be solved in terms of Jacobian elliptic 
functions by applying a substitution which converts the problem to the form 
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(^) = (1 ( 1 ”^) 

This substitution depends on the degree of f*? 4 and its number of complex zeros as 
well as on the type of orbit you want to obtain. For a list see e.g. IB. 

To find a solution of (fT2l i in terms of the Weierstrass elliptic function you have 
to convert the problem to the standard form (flSl l. which can be obtained by first 
converting a polynomial of degree four to degree three by x = +4C/> if f ’3 4 = 

and F 3 4 (x/>) = 0 and subsequently, or if P 3 4 was of degree three in the first 
place, substituting ^ = ^(4z+ y) if P = is the polynomial of degree three. 

Note that you always have to take care of the initial condition, too. 


3.2 Solutions in terms of hyperelliptic functions 


To generalize the solution methods outlined in the previous section we first need to 
consider the Jacobi inversion problem 


yi = 



t' ^dt 

TWV 


/=i,...,g 


( 18 ) 


where P{t) = is a polynomial of degree 2 .g + 1 and g is the genus. 

Note that forg = 1 and ^2 = 0 we recover (ffSl l. The g solutions X;(yi,... ,yg) of (fT^ 
can be given in terms of generalized Weierstrass functions. These are defined by the 
theta function via the Kleinian sigma function 

(19) 

pij{z-, ®i, ® 2 ) = - ^ ^ log a(z; ®i, ® 2 ), ( 20 ) 

OZi OZj 

where z € C^, O); are g x g matrices such that Wf ^©2 is symmetric with positive 
definite imaginary part, K = 77 ( 2 ©i)^* with the periods of the second kind 2rj, 
Koc is the vector of Riemann constants, and C is a constant which does not matter 
here (for further details see e.g. 121; note thatK<x. = 5 )' + ( 0 , 5 )^ if g = 2). The 

solutions of ( fTSb are then given by the solutions of 


+ J^Pgi{yi,---,yg)x‘ ^ = 0 . 


!=1 


( 21 ) 


Let us consider now the case g = 2 and a general differential equation of the form 


dx 

dy 


^5.6 W, ■*(yo)=-*o 


( 22 ) 
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where P 5 5 is a polynomial of degree five or six. If it is of degree six it can be 
reformulated as {x) by a substitution x = +xp, where xp is a zero 

of P 5 6 and is a polynomial of degree five. Such a differential equation, or, if P 5 g 
was of degree five in the first place, can then be cast in the form 

= i=h2 (23) 

ay 

by an appropriate normalization. We may then find the solution to this equation as 
the limiting case X 2 —of the Jacobi inversion problem (fTSl l in the following way: 
first observe that 


t :=x\= lim 


= lim 


X\X2 
+X2 

(7(Ji2 — <7iCJ2 
(7,2 - (7(722 


= lim 


Pniyuyi) 


piiiyuyi) 

(yuyi), 


(24) 


where ( 7 ,- (z) is the derivative of (7 with respect to z,-. From (fTSl l with g = 2 and X 2 —?-°° 
we may identify either yi or yz with our physical coordinate y in the differential 
equation (|2^ . say y,. Fortunately, we get automatically rid of the other yj, j ^ i, 
by the same limiting process X 2 This is because the set of zeros of the theta 
function z 1 —0 [/Cx,]((2(b)^^z), which is a one dimensional submanifold of C^, is 
given by all vectors z = (zi , 22 ) which can be written as z, = 


•V t'-'dt 


with the same 


X (see e.g. ED). This is exactly true for the vector (yi,y 2 ), which means that we 
may write yj = f{yi) for some function /. As the zeros of the theta function are also 
zeros of (7 we can simplify (l24l i to 


t = -^iy\J{y\)) 

C72 


t = -^{f(y2),y2)- 

02 


(25) 


Note that according to the given initial condition we actually have that y,- is the 
physical coordinate minus a constant. 


4 Analytical solutions in the literature 

In this section we will collect applications of the methods outlined in section [3 to 
geodesic motion in the Kerr-Newman-Taub-NUT-de Sitter spacetime as given in 
section 12] For older literature we refer to Sharp llTTl who collected most of the pa¬ 
pers on geodesic motion in Kerr-Newman spacetime and subclasses, which were 
available at that time. Partly this is still quite complete but we also try to update his 
collection (with respect to analytical solutions). Note that we only consider analyti¬ 
cal solutions to general timelike and lightlike geodesics (with an electric or magnetic 
charge, as applicable). In particular, we do not list the vast literature on equatorial 
motion in Kerr spacetime. Of course, we do not claim that our list is complete. 
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Schwarzschild: Regarding analytical solution methods the list of Sharp already 
contained the complete set of solutions. Most notably, this includes the works by 
Hagihara na, who derived the analytical solutions in terms of Weierstrass elliptic 
functions, and Darwin iiiiia, who used Jacobian elliptic functions and integrals. 

Reissner-Nordstrom: Surprisingly, the analytical solutions to the geodesic equa¬ 
tion in Reissner-Nordstrom spacetime seem to be considered first only in 1983 by 
Gackstatter |l2l although it can be handled completely analogously to the Schwarzschild 
case. He studied bound timelike geodesics and light in terms of Jacobian elliptic in¬ 
tegrals and functions. Recently, Slezakova ll28l gave a comprehensive analysis of ar¬ 
bitrary timelike, lightlike, and even spacelike geodesics. Grunau and Kagramanova 
m solved the equations of motion of electrically and magnetically charged parti¬ 
cles in Reissner-Nordstrom spacetime in terms of Weierstrass elliptic functions. 

Taub-NUT: Timelike geodesics were studied by Kagramanova et al in terms 

of Weierstrass elliptic functions. 

Kerr: Most of the older literature on Kerr spacetime is concerned with the much 
simpler particular case of equatorial geodesics. We refer to Sharp ET) here for these 
early works. Note that in terms of the proper time (or the corresponding affine pa¬ 
rameter for light) the equations of motion are still coupled. Therefore, most of the 
analytical solutions before the introduction of the Mino time ll25]l implicitly included 
integrals over the latitude or the radius, see e.g. Kraniotis or Slezakova ll28l for 
a review. As notable exception, Cadez et al 1291 introduced already in 1998 a similar 
parameter (called P, see their equation (34)) as they considered the motion of light. 

After the introduction of the Mino time, in 2009 Fujita and Hikida @ used this 
new affine parameter to derive the analytical solution for bound timelike geodesics 
in terms of Jacobian elliptic functions. General timelike geodesics and lightlike mo¬ 
tion were treated shortly after that by Hackmann HU in 2010. 

Note that Kraniotis l22l also derived analytical solutions for lightlike geodesics 
in terms of hypergeometric functions. 

Kerr-Newman: Charged particle motion was considered by Hackmann and Xu 
ini in terms of Weierstrass elliptic functions. 

Schwarzschild-de Sitter, also called Kottler space-time; Note that on the level of 
the differential equation, lightlike geodesics in Schwarzschild-de Sitter are identical 
with the lightlike equations of motion for Schwarzschild, as the cosmological con¬ 
stant can be absorbed in the definition of just a single parameter. Analytical solution 
are given e.g. in Gibbons et al HI. General timelike geodesics in Kottler spacetime 
can be treated in terms of hyperelliptic functions as elaborated by Hackmann and 
Lammerzahl nma. 

Reissner-Nordstrom-de Sitter: The equations of motion for general timelike 
geodesics were solved in ca. The motion of photons was very recently analyti¬ 
cally calculated by Villanueva et al for a negative cosmological constants using 
Weierstrass elliptic functions. 

Taub-NUT-de Sitter: Time like motion was analyzed in lfT3l . 

Kerr-de Sitter: The equations of motions for timelike geodesics were analyti¬ 
cally solved by Hackmann et al uniiibi in terms of hyperelliptic functions. Note 
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that Kraniotis 1221 also derived analytical solutions for lightlike geodesics in terms 
of hypergeometric functions. 

Kerr-Newman-Taub-NUT-de Sitter: The general solution for timelike geodesic 
was shortly outlined in 1131 . 

Acknowledgements We thank the German Research Foundation DFG for financial support within 
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References 


1. M. Ahramowitz and I. A. Stegun, editors. Handbook of Mathematical Functions. Dover 
Publications, Inc., New York. 

2. V.M. Buchstaber, V.Z. Enolskii, and D.V. Leykin. Hyperelliptic Kleinian Functions and Ap¬ 
plications. Reviews in Mathematics and Mathematical Physics 10. Gordon and Breach, 1997. 

3. B. Carter. Phys. Rev., 174:1559, 1968. 

4. C. Darwin. Proc. R. Soc. London A, 249:180, 1959. 

5. C. Darwin. Proc. R. Soc. London A, 263:39, 1961. 

6. R. Eujita and W. Hikida. Class. Quantum Grav, 26:135002, 2009. 

7. E. Gackstatter. Ann. Phys., 495:352, 1983. 

8. G. W. Gibbons, C. M. Wamick, and M. C. Werner. Class. Quantum Grav., 25:245009, 2008. 

9. J.B. Griffiths and J. Podolsky. Int. J. Mod. Phys., 15:335, 2006. 

10. S. Gmnau and V. Kagramanova. Phys. Rev. D, 83:044009, 2011. 

11. E. Hackmann. Geodesic equations in black hole .space-times with cosmological constant. PhD 
thesis, Universitat Bremen, 2010. 

12. E. Hackmann, V. Kagramanova, J. Kunz, and C. Lammerzahl. Phys. Rev. D, 78:124018, 2008. 

13. E. Hackmann, V. Kagramanova, J. Kunz, and C. Lammerzahl. Europhys. Lett., 88:30008, 
2009. 

14. E. Hackmann and C. Lammerzahl. Phys. Rev. Lett., 100:171101, 2008. 

15. E. Hackmann and C. Lammerzahl. Phys. Rev., D 78:024035, 2008. 

16. E. Hackmann, C. Lammerzahl, V. Kagramanova, and J. Kunz. Phys. Rev, D 81:044020, 2010. 

17. E. Hackmann and H. Xu. Phys. Rev. D, 87:124030, 2013. 

18. Y Hagihara. Japanese J. Astr. Geophys., 8:67, 1931. 

19. Adolf Hurwitz. Vorlesungen tiber Allgemeine Funktionentheorie und elliptische Funktionen. 
Springer-Verlag, Berlin, 1964. 

20. V. Kagramanova, J. Kunz, E. Hackmann, and C. Lammerzahl. Phys. Rev., D 81:124044, 2010. 

21. G.V. Kraniotis. Class. Quantum Grav, 22:4391, 2005. 

22. G.V. Kraniotis. Class. Quantum Grav, 28:085021, 2011. 

23. V. S. Manko and E. Ruiz. Class. Quantum Grav, 22:3555, 2005. 

24. A. 1. Markushevich. Theory of functions of a complex variable. Vol. I, If HI. Chelsea Pub¬ 
lishing Co., New York, english edition, 1977. Translated and edited by Richard A. Silverman. 

25. Y. Mino. Phys. Rev, D 67:084027, 2003. 

26. D. Mumford. Tata Lectures on Theta, Vol. I and IT Birkhauser, Boston, 1983/84. 

27. N. A. Shaip. Gen. Relativ. Gravitation, 10:659, 1979. 

28. G. Slezakova. Geodesic Geometry of Black Holes. PhD thesis, University of Waikato, 2006. 

29. A. Cadez, C. Eanton, and M. Calvani. New Astronomy, 3:647, 1998. 

30. J.R. Villanueva, J. Saavedra, M. Olivares, and N. Cruz. Astrophys. Space Sci., 344:437, 2013. 


